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L. Ya. GlozmanQ C. B. LangQ and M. Schr6cl|] 

Institut fur Physik, FB Theoretische Physik, Umversitat Graz, A-8010 Graz, Austria 

(Dated: March 4, 2013) 

We study hadron correlators upon artificial restoration of the spontaneously broken chiral symme- 
try. In a dynamical lattice simulation we remove the lowest lying eigenmodes of the Dirac operator 
from the valence quark propagators and study evolution of the hadron masses obtained. All mesons 
and baryons in our study, except for a pion, survive unbreaking the chiral symmetry and their 
exponential decay signals become essentially better. From the analysis of the observed spectro- 
scopic patterns we conclude that confinement still persists while the chiral symmetry is restored. 
All hadrons fall into different chiral multiplets. The broken U(1)a symmetry does not get restored 
upon unbreaking the chiral symmetry. We also observe signals of some higher symmetry that in- 
cludes chiral symmetry as a subgroup. Finally, from comparison of the A — TV splitting before 
and after unbreaking of the chiral symmetry we conclude that both the color-magnetic and the 
flavor-spin quark-quark interactions are of equal importance. 
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I. INTRODUCTION 

Highly excited hadrons in the u, d sector reveal some 
parity doubling [TUTU] and possibly some higher symme- 
try. It was conjectured that this parity doubling reflects 
effective restoration of chiral symmetry, i.e., insensitivity 
of the hadron mass generation mechanism to the effects of 
chiral symmetry breaking in the vacuum 1-6 . Whether 
this conjecture is correct or not can be answered exper- 
imentally since the conjectured symmetry requires exis- 
tence of some not yet observed states. 

Recent and most complete experimental analysis on 
highly excited nucleons that includes not only elastic ttN, 
but also the photoproduction data, does report evidence 
for some of the missing states and the parity doubling 
patterns look now even better than before |llj . 

The question of a possible symmetry in hadron spectra 
is one of the central questions for QCD since it would 
help to understand dynamics of confinement and chiral 
symmetry breaking as well as their role for the hadron 
mass generation. 

Another "experimental" tool to address the issue of the 
hadron mass generation is lattice QCD. Equipped with 
the QCD Lagrangian and Monte-Carlo techniques, one 
can calculate, at least in principle, hadron masses and 
other hadron properties from first principles. Enormous 
progress has been achieved for the hadron ground states. 
The problem of excited states, especially above the mul- 
tihadron thresholds like irN, Ait, 7T7t, irp, . . . turns out 
to be much more difficult and demanding than was ini- 
tially anticipated. When it is solved lattice results should 
reproduce experimental patterns and possibly indicate 
some still missing states. 
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Still, the mass of a hadron by itself, obtained from the 
experiment or from the lattice simulations, tells us not so 
much about the physics which is behind the mass gener- 
ation. The pattern of all hadrons, on the contrary, could 
shed some light on the underlying dynamics if there are 
some obvious symmetries in the pattern or if its regular- 
ities can be systematically explained. 

The most interesting issue is to get some insight on 
how QCD "works" in some important cases and under- 
stand the underlying physical picture. In this sense one 
can use lattice QCD as a tool to explore the interrela- 
tions between confinement and chiral symmetry break- 
ing. In particular, we can ask the question whether 
hadrons and confinement will survive after having artifi- 
cially removed the quark condensate of the vacuum. This 
can be achieved via removal of the low-lying eigenmodes 
of the Dirac operator, which is a well defined procedure 

[una. 

In the past mainly the opposite was explored. After 
suggestions within the instanton liquid model |14) the 
effect of the low-lying chiral modes on the p and other 
correlators was studied on the lattice. In a series of pa- 
pers O [T5l - [17] it was shown that low modes saturate 
the pseudoscalar and axial vector correlators at large dis- 
tances and do not affect the part where high-lying states 
appear. In [TH [TH] low mode saturation and also ef- 
fects of low mode removal for mesons were studied for 
quenched configurations with the overlap Dirac opera- 
tor |20] • Subsequently low modes were utilized to 
improve the convergence of the determination of hadron 
propagators (TSJ, HH [2TH24] studying the efficiency when 
using the low modes of the Dirac operator or the Hermi- 
tian Dirac operator. 

We are studying the complementary case, i.e., removal 
of the low modes and we will refer to this as "unbreak- 
ing" the chiral symmetry. This issue has been addressed 
in a recent paper 1251 126] where the low- lying eigenmodes 
of the Dirac operator have been removed from the quark 
Green's function and masses of the lowest mesons tt, p, ap 
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and ai have been calculated with such truncated quark 
propagators. The truncated Landau gauge quark propa- 
gator itself has been investigated in [27J where the loss of 
dynamical mass generation in the infrared sector of the 
propagator has been demonstrated. 

After the unbreaking of the chiral symmetry the signal 
from the 7r-meson, obtained with the pseudoscalar quark- 
antiquark operator, disappears, which is consistent with 
the (pseudo) Goldstone boson nature of the pion. In- 
deed, with the artificially restored chiral symmetry there 
cannot be Goldstone bosons. What is very interesting, 
is that other low-lying mesons survive and the quality of 
their signals become even essentially better after extrac- 
tion of the low-lying eigenmodes of the Dirac operator, 
responsible for chiral symmetry breaking. The very fact 
that hadrons survive the unbreaking of the chiral sym- 
metry tells that there is confinement in the system even 
without the quark condensate. (A similar behavior was 
found in [28j ESI , where the effect of such a removal on 
the static quark potential was studied.) 

After extraction of the quark condensate the lowest- 
lying p and a\ mesons demonstrate restoration of the 
chiral symmetry - they become degenerate - and their 
mass is rather large. This disproves a rather popular as- 
sertion that, e.g., the p-mass is entirely due to the quark 
condensate of the vacuum. As a physical implication we 
should then not expect a drop off the p-mass in a dense 
medium, which is a very popular issue both theoretically 
and experimentally for the last two decades [3D]. This 
result should also be of importance for a debated issue of 
confining but chirally symmetric matter at low tempera- 
tures and large density [3TH34] . 

The conclusion of |25j about survival of confinement 
after unbreaking of the chiral symmetry was obtained 
on the limited basis of the lowest-lying mesons. In or- 
der to see it more clearly we need to extend the number 
of extracted states, in particular to include radially and 
orbitally excited hadrons (following quark model termi- 
nology). Consequently, we now add the b\ and p' states 
to the above mentioned list of mesons. We limit ourselves 
only to the isovector mesons since the isoscalar mesons 
would require inclusion of disconnected graphs, which is 
numerically very costly. 

Most importantly, we study behavior of the ground 
and excited positive and negative parity states in the N 
and A spectra. The baryonic states add additional infor- 
mation about existence or nonexistence of confinement. 
They do allow to see that confinement does survive af- 
ter the restoration of chiral symmetry. Second, we see 
some traces of the higher symmetry, higher than simply 
SU(2)r x SU(2)l. This observation may be related with 
the higher symmetry seen in the highly excited hadrons. 

Our paper is organized as follows. In the next section 
we remind on the connection between low eigenmodes of 
the Dirac operator and the quark condensate in the vac- 
uum. We discuss some basic aspects of removal of Dirac 
eigenmodes from the quark propagator. In Sect. |III| we 
present the details of the lattice simulation. The fourth 



section is devoted to the description of the baryon and 
meson interpolators used for our study of hadrons. In 
Sect. [V]we show and analyze the results and draw con- 
clusions. Finally, in the last section, we briefly summarize 
our main observations. 



II. THE QUARK CONDENSATE AND THE 
DIRAC OPERATOR 

The lowest eigenmodes of the Dirac operator are re- 
lated (in the chiral limit) to the quark condensate of the 
vacuum. This is encoded in the Banks-Casher relation 

m 

< 0\qq\0 >= -7T P (0) , (1) 

where p(0) is a density of the lowest quasi-zero eigen- 
modes of the Dirac operator. Here the sequence of lim- 
its is important: first, the infinite volume limit at finite 
quark mass is assumed and then the chiral limit should 
be taken. The opposite sequence would produce no chiral 
condensate as in the finite volume there cannot be any 
spontaneous breaking of chiral symmetry. 

All lattice calculations are performed on a lattice of 
a finite volume. In the finite lattice volume the spec- 
trum of the Dirac operator is discrete and the energy of 
the lowest nonzero mode of the Dirac operator is finite. 
Consequently, the quark condensate is strictly speaking 
zero. However, increasing the lattice volume the gap 
in the spectrum of the eigenmodes becomes smaller and 
smaller and the density of the lowest nonzero eigenmodes 
increases. A well defined limit of this density scaling ex- 
ists: the number of such eigenvalues in a given interval 
adjacent to the real axis scales with the lattice volume. 
(In |36| is is arguecQ that the number of relevant eigen- 
values should scale proportional the square root of the 
number of lattice points.) 

In [24j[25] it was established that the eigenmodes of the 
Hermitian Dirac operator D 5 = 75 D result in a faster 
saturation of the pseudoscalar correlator when approxi- 
mating quark propagators by the lowest eigenmodes only, 
compared to the eigenmodes of the Dirac operator D. 
Therefore, we focus on reducing the quark propagators 
in terms of eigenmodes of D5 rather than D. 

From the lattice calculations in a given finite volume 
we cannot say a priori which and how many lowest eigen- 
modes of the Dirac operator are responsible for the quark 
condensate of the vacuum. For the overlap operator the 
real eigenvalues correspond to exact chiral modes, the 
zero modes (instantons); for Wilson- type operators one 
may associate real eigenvalues with zero modes. Their 
weight is suppressed in the infinite volume limit. The 
Banks-Casher relation, however, relies only on the den- 
sity of nonzero modes. For the Hermitian Dirac operator 
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D5 there is no simple method to distinguish between the 
real modes of the D and its small, but complex modes. 
In [251 US] we discuss the integral over the distribution 
of the (real) eigenvalues of D$. There we observe a tran- 
sition region up to roughly twice the size of the quark 
mass corresponding to 0(16 — 32) eigenmodes, as also 
observed in, e.g., [3"7H4"0] . 

We follow the procedure to remove an increasing num- 
ber of the lowest Dirac modes and study the effects of 
the (remaining) chiral symmetry breaking on the masses 
of hadrons. To be specific, we construct reduced quark 
propagators 

S Ie d(k) = S - Si m ( fe ) = S -^2^ 1 \v l )(v l \-f 5 , (2) 

i<k 

where S is the standard quark propagator obtained from 
the inversion of the Dirac operator, the pn are the (real) 
eigenvalues of D5, \vi) are the corresponding eigenvectors 
and k represents the reduction parameter which will be 
varied from — 128. 

Note that the low-mode contribution (lm) of the quark 
propagators, Swfc), must act on the same quark sources 
as S, see discussion below. 

III. THE SETUP 
A. Dirac operator 

For the dynamical quarks of our configurations as well 
as for the valence quarks of our study the so called chi- 
rally improved (CI) Dirac operator [311 [32] has been 
used. The latter represents an approximate solution to 
the Ginsparg- Wilson equation and therefore offers better 
chiral properties than the Wilson Dirac operator while 
being less expensive, in terms of computation time, in 
comparison to the chirally exact overlap operator. 

B. Gauge configurations 

We performed our study on 161 gauge field configu- 
rations [321 [33] that were generated for two degenerate 
dynamical light CI fermions with a corresponding pion 
mass m v = 322(5) MeV. The lattice size is 16 3 x 32 and 
the lattice spacing a = 0.144(1) fm. 

C. Quark source smearing 

In order to obtain quark propagators, the Dirac oper- 
ator has to be inverted on given quark sources. To im- 
prove the signal in hadron correlators, extended sources 
of Gaussian form (45j 06] instead of point sources are 
being used. Using several different extended sources al- 
lows for a larger operator basis in the variational method 
[3ZH3S]- We use three different kinds of sources: narrow 



(0.27 fm) and wide (0.55 fm) sources, which are approxi- 
mately of Gaussian shape, and a derivative source. 

The narrow (wide) sources will be denoted by a sub- 
script n (w) of the quark fields and the derivative source 
by di, respectively. The details of the calculation of the 
smeared quark sources are given in |44j . 

D. Variational method 

In order to disentangle the excited states from the 
ground state (and also to provide cleaner signals for the 
ground states) we use the variational method [33 EH]- 
One computes cross-correlators Ctk{t) = (Oi(t)Ok(0y) 
between several different lattice interpolators and solves 
the generalized eigenvalue problem 

C(t)u n (t) = X n (t)C(t )Un(t) , (3) 

in order to approximately recover the energy eigenstates 
|n). The eigenvalues allow us to get the energy values 
X n (t) exp(— E n t) and the eigenvectors serve as finger- 
prints of the states, indicating their content in terms of 
the lattice interpolators. In our plots we show X n (t), the 
effective masses E n (t) = log(A n (t)/A„(t + 1)) and the t- 
dependence of the eigenvectors in order to verify the state 
identification. The quality of the results depends on the 
statistics and the provided set of lattice operators (for a 
discussion see [44]). The used interpolators are discussed 
in Sect. El 

E. Dirac eigenmodes 

On the given gauge field configurations we calculated 
the lowest 128 eigenmodes of the Hermitian Dirac oper- 
ator L>5 using ARPACK which is an implementation of 
the Arnoldi method to calculate a part of the spectrum 
of arbitrary matrices [50] . 

Once the eigenmodes have been calculated and the 
quark propagators S have been obtained by inverting the 
Dirac operator on the three types of sources mentioned 
in Sec. |III C[ we can construct the reduced propagators 
•Sred(fe) according to (|2j) after multiplying the low-mode 
part of the propagator, S\ m tk) , with the same three source 
types, respectively. 

IV. HADRON INTERPOLATORS 

Here we list the baryons and mesons we studied un- 
der Dirac low-mode reduction and give the interpolating 
fields for each individual. 

A. Baryons 

We analyze the nucleon and A baryons both with pos- 
itive and negative parity. For the interpolators we use 
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TABLE I. Interpolators for the N channel. The Dirac struc- 
tures % j the quark smearings and the corresponding inter- 
polator numbers #jv are given. 



Gaussian smeared quark sources (n and For the nu- 
cleon we adopt three different Dirac structures, resulting 
in f 8 interpolators (see Tab. [I]) where we left out those 
operators that are similar to other ones due to isospin 
symmetry. The construction of the nucleon interpolators 
is given by 



=e„ 6c lf ) u a (ulTf d, 



(i) 



For the A, 



A fc = e abc u a (uj C 7 fc u c ) 



(4) 



(5) 



we use only one Dirac structure and the six corresponding 
interpolators are listed in Tab. [TT] We use parity projec- 
tion for all baryons and Rarita-Schwinger projection for 
the A [33] . The sink interpolators are also projected to 
zero spatial momentum. 
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TABLE II. Interpolators for the A channel. The quark smear- 
ings and the corresponding interpolator numbers #a are 
given. 
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TABLE III. Interpolators for (top) the p-meson, J 1 



(middle) the ai-meson, J 
meson, J PC — 1 



= 1 , and (bottom) the b\- 
The first column shows the number, 
the second shows the explicit form of the interpolator. The 
numbers refer to the classification in 1521. 



operator basis for the variational method. In Table III 
we list only those interpolators explicitly whose combi- 
nation resulted in a good signal in practice when plugged 
into the variational method. The sink interpolators are 
projected to zero spatial momentum. A more complete 
list of possible interpolating fields is given in [33J ■ 



V. RESULTS AND DISCUSSION 



B. Mesons 

We investigate isovector mesons of spin 1. Isoscalars 
require the calculation of disconnected graphs which are 
computationally too demanding for the type of fermion 
action used. The scalar meson ao as well as the pseu- 
doscalar pion was studied already in [25 . 

Thus, the studied nonexotic channels are the J com- 
binations 1 (p), 1 ++ (ai) and l" 1 (&i). For the ana- 
lysis of the mesons we include derivative sources |51| in 
the construction of the interpolators to provide a large 



A. Truncation study 

The reduction parameter k in ([2| gives the number 
of the lowest eigenmodes of the Dirac operator removed 
from the quark propagator. We study reduced quark 
propagators with k = 0,2,4,8,12,16,20,32,64,128 and 
three different quark smearings (n, w, di). These are 
then combined into different hadron propagators and the 
correlation matrix for a hadron with given quantum num- 
bers is calculated. The variational method is used to ex- 
tract the ground and excited states of that hadron. Con- 
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FIG. 1. p with 12 eigenmodes subtracted: The correlators for all eigenstates (upper left), effective mass plot for the two lowest 
states (upper right), eigenvectors corresponding to the ground state (lower left), and 1st excited state (lower right). 
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FIG. 2. p with 32 eigenmodes subtracted: The correlators for all eigenstates (upper left), effective mass plot for the two lowest 
states (upper right), eigenvectors corresponding to the ground state (lower left), and 1st excited state (lower right). 



sequently, we observe and study the evolution of hadron 
masses as a function of the number of the subtracted low- 
est eigenmodes. Increasing the number of the subtracted 
lowest eigenmodes we gradually remove the chiral con- 
densate of the vacuum and consequently "unbreak" the 
chiral symmetry. 

Typical results for the mesons and baryons under study 
are shown in Figs. [T] - [l2j For each hadron we show 
in the figures two representative reduction levels k for 
which we also show explicitly all eigenvalues stemming 
from the variational method, i.e., the correlators corre- 
sponding to different energy levels. Moreover, we show 
for the ground and first excited state (where applicable) 
the eigenvector components and the effective mass plots 
including fit ranges and values. The energy values are 
determined from exponential fits to the eigenvalues over 
the indicated fit ranges. 

In Figs. [T] a nd [2] w e show the eigenvalues and eigenvec- 
tors (cf. Sect. |III D ) and the effective mass plots for the 



ground and excited states of the p- meson (J PC = 1 ) 
after having subtracted 12 and 32 eigenmodes of D5. 
(The results for the untruncated situation, with higher 
statistics and for several more parameter sets are shown 
in |52j.) The eigenvector composition for both states is 
stable and clearly distinct. The mass splitting between 
ground state and excited state disappears with increasing 
truncation level. 

For the meson channels a\ and b\ the statistics allow 
only to determine the ground state in a reliable way. In 
the case of the a\ (J PC = 1 ++ ) meson (Fig. [3] for k — 4 
and k = 64) and for the bi {J PC = l+~) meson (Fig. [4] 
for k — 2 and k = 128) we thus show only the ground 
states. We observe improved plateau quality for the effec- 
tive masses when increasing number of truncated modes. 

In Fig. [5] the nucleon of positive parity (the ground 
state and its first excitation, J p = ^ + ), after having 
subtracted the lowest 20 eigenmodes is shown. These 
should be compared with Fig. [6] where 64 modes have 
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FIG. 4. 61 with 2 (upper row) and 128 (lower row) eigenmodes subtracted: The correlators for all eigenstates (left), effective 
mass plot for the lowest state (right). 



been subtracted. In Fig. [7] and [8] we present two nucleon 
states of negative parity, J P — \ , (at reduction level 12 
and 64). 

The positive parity A ground and excited states, J p = 
| + at reduction level 16 and 128 are shown in Fig. [9] and 
Fig." 
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and the negative parity A's (J p = | ) at the 
same reduction levels are given in Fig. [TT] and Fig. [T2| 
respectively. 

An obvious observation is that the quality of the sig- 
nal (the quality of plateaus) essentially improves with 
increasing the number of removed eigenmodes for all 
hadrons under study. This fact makes it easier to reli- 
ably identify masses of states after unbreaking the chiral 
symmetry. In many cases for excited hadrons the qual- 
ity of plateaus before unbreaking of the chiral symmetry 
is rather poor, but after removing more and more eigen- 
modes of the Dirac operator it becomes better and better, 
so eventually it allows to unambiguously establish that 
indeed we see the state, even though with the untrun- 



cated propagators the identification of the state would 
be less clear. 

A plausible explanation for this phenomenon would be 
that by unbreaking the chiral symmetry we remove from 
the hadron its pion cloud and subtract all higher Fock 
components like 7riV, 7rA, tttt, and so on from the hadron 
wave functions. It is these components related to the 
chiral symmetry breaking that couple excited and not 
excited hadrons to each other and render signals from 
the excited states poor in fully untruncated QCD. 



B. Confinement after unbreaking the chiral 
symmetry 

The most interesting question is whether hadrons and 
confinement survive the unbreaking of the chiral sym- 
metry. To discuss this issue we put all the results from 
the previous figures together and analyze how masses of 
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the considered hadrons change with the reduction level. 
Therefore, the relevant truncation scale is not the Dirac 
operator eigenmode index k itself, since it has to scale 
with the lattice volume when keeping the physics con- 
stant. Instead, we introduce a cutoff parameter a such 
that the reduction level a means that all jik for which 
\fj,f.\ < a have been excluded in the underlying quark 
propagators [25j . We still give the corresponding index k 
on the upper horizontal scale of the plots. 

The masses of all studied hadrons under D 5 eigenmode 
reduction are summarized in Fig. |13| The scale is set by 
the Sommer parameter (i.e., by the static potential act- 
ing between two heavy quarks, in other words, by the 
gluonic dynamics). The physics of the static potential 
knows nothing about the valence quarks and the trun- 
cation of lowest lying eigenmodes of the Dirac operator 
for the light valence quarks. It implies that these plots 
suggest at least qualitatively the evolution of the hadron 
masses in absolute units of energy. We observe approxi- 
mately a universal growth of all hadron masses with equal 
slope after subtraction of a sufficient amount of the chiral 
modes of the Dirac operator. In this regime chiral sym- 
metry is approximately restored and the whole hadron 
mass is not related to chiral symmetry breaking. A uni- 
versal slope might be interpreted as an indication to a 
universal growth of the hadron size. 

Assume that after having unbroken the chiral symme- 
try the exponential decay signals from all hadrons would 
disappear. This would indicate that with the artificial 
restoration of the chiral symmetry confinement also van- 
ishes and that there is a direct connection between the 
confinement in QCD and the lowest lying modes of the 
Dirac operator. Contrary to that we observe a very clear 
signal from all hadrons, except for a pion. This suggests 
that confinement survives the unbrcaking of the chiral 
symmetry. 

However, there is still the possibility that this clear 
signal comes from the unconfined (unbound) quarks with 
some mass mo at a given truncation level. In this case 
we would expect at a given truncation level a universal 
scaling law with all mesons having the mass 2 mo and all 
baryons with the mass 3m . There would be no excited 
states of hadrons. 

In order to address this issue we show in Fig. [14] all 
hadron masses in units of the p-meson mass obtained at 
the same truncation level. Indeed, some of the states - 
such as ai and (/, as well as the ground and the first 
excited states of the nucleon of both parities - do fol- 
low this behavior of mass 2 mo and 3 m for mesons and 
baryons, respectively. However, this is definitely not the 
case for the b\ state as well as for the A-resonance and 
especially its first excited states of positive and negative 
parity. Given that the signal in all latter cases is unam- 
biguous, we conclude that there is no universal scaling 
(2 mo for mesons and 3 mo for baryons) for all hadrons. 
This rules out the possibility that our signals are pro- 
duced by the unbound (unconfined) quarks. We do ob- 
serve confined hadrons. 



Actually, this can be seen also from another perspec- 
tive. The mass mo is large and increases with truncation 
of the quark propagators. At the same time we observe 
chiral restoration in the correlators (e.g., a\ and p). The 
large mass mo of unconfined free quarks then contradicts 
restoration of chiral symmetry. This supports our argu- 
ment that we do not observe unconfined quarks. 

The fact that masses of some of the mesons and some 
of the baryons get degenerate and are related through 
a simple law 2 mo for mesons and 3 mo for baryons indi- 
cates symmetries of hadrons that appear after unbreaking 
(restoration) of the chiral symmetry. 



C. Meson degeneracies and splittings and what 
they tell us 

Restoration of the SU(2)l x SU(2)r chiral symmetry 
in the vacuum requires the states to fall into parity-chiral 
multiplets [2HS]- Below we shortly summarize a content 
of these parity-chiral multiplets for the mesons in our 
study. 

The full set of multiplets of the parity-chiral group 
SU(2) L x SU(2) R x d, where the group Cj consists of 
identity and the space inversion, for the J = 1 mesons is 
as follows: 



(0,0) 

(|j §)a 
(si k)b 

(0,1) + (1,0) 



w(0,l— ) 
fci(0,l + ") 

w(o,i—) 

oi(l,l++) 



A(o,i++) 
) 

Oi(U + -) 



Note, that the unbroken chiral SU(2)l x SU(2)r sym- 
metry requires existence of two independent p-mesons, 
one of them is the chiral partner of the h\ meson, and 
the other one of the ai state. Similar is true for the 
w-meson. 

The states from two distinct multiplets (|, \ ) a and 
(|, i)f, that have the same isospin but opposite spatial 
parity are connected to each other by the U{1)a trans- 
formation, if the U(1)a symmetry is broken neither ex- 
plicitly nor spontaneously. In our real world U(1)a is 
broken both explicitly via the axial anomaly and sponta- 
neously via the quark condensate of the vacuum. So in 
the world with restored U(1)a symmetry a p meson, that 
is the chiral partner to the hi meson, would be degenerate 
with the b\ state. The hi, p, w and b\ states would form 
an irreducible multiplet of the SU{2) L x SU(2) R x U(1) A 
group. 

On top of the chirally symmetric vacuum the p — a\ 
splitting vanishes, see, e.g., Figs. [T~3| and [H} , a clear sig- 
nal of the chiral SXJ{2)l x SU(2)p symmetry restoration 
in the physical states. At the same time large b\ — p 
and b\ — p' splittings persist. This is a direct indication 
that the U(1)a breaking does not disappear. While that 
U(1)a breaking component that is due to the chiral con- 
densate should vanish with the condensate, the U(1)a 
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FIG. 5. N(+) with 20 eigenmodes subtracted: The correlators for all eigenstates (upper left), effective mass plot for the two 
lowest states (upper right), eigenvectors corresponding to the ground state (lower left), and 1st excited state (lower right). 
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FIG. 8. N(-) with 64 eigenmodes subtracted: The correlators for all eigenstates (upper left), effective mass plot for the two 
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FIG. 9. A(+) with 16 eigenmodes subtracted: The correlators for all eigenstates (upper left), effective mass plot for the two 
lowest states (upper right), eigenvectors corresponding to the ground state (lower left), and 1st excited state (lower right). 





FIG. 10. A(+) with 128 eigenmodes subtracted: The correlators for all eigenstates (upper left), effective mass plot for the two 
lowest states (upper right), eigenvectors corresponding to the ground state (lower left), and 1st excited state (lower right). 
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FIG. 11. A(— ) with 16 eigenmodes subtracted: The correlators for all eigenstates (upper left), effective mass plot for the two 
lowest states (upper right), eigenvectors corresponding to the ground state (lower left), and 1st excited state (lower right). 
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FIG. 12. A(— ) with 128 eigenmodes subtracted: The correlators for all eigenstates (upper left), effective mass plot for the two 
lowest states (upper right), eigenvectors corresponding to the ground state (lower left), and 1st excited state (lower right). 



breaking via the axial anomaly still persists. Then it fol- 
lows that there is no direct interconnection of the lowest 
lying modes of the Dirac operator and the mechanism 
of the anomalous U{1)a breaking in QCD. Such a direct 
interconnection was suggested in the past through, e.g., 
the instanton fluctuations. 

After unbreaking of the chiral symmetry the p and p' 
mesons become degenerate. What does this tell us? A 
degeneracy indicates some symmetry. The two distinct p 
states, p and p' , lie in different irreducible parity-chiral 
representations, (|, \) a and (0,1) + (1,0). In principle, 
their degeneracy could point out to a reducible represen- 
tation of the parity-chiral group that would include both 
irreducible representations. Indeed, the product of two 
fundamental quark-antiquark chiral representations does 
contain, in particular, both (±, i) a and (0, 1) + (1,0): 



Such a multiplet of dimension 16 (including isospin de- 
generacies) would consist of two distinct uj- mesons, fx, 
hi, two p-mesons as well as b\ and a\ mesons and would 
require a degeneracy of all of them. Now we do find, 
however, that the b\ meson is well split from both p 
and p' after the unbreaking of the chiral symmetry. This 
rules out that the observed p — p' degeneracy is related 
to restored chiral symmetry. The degenerate p and p' 
states are different because their eigenvectors are orthog- 
onal and because they are well split before the removal 
of the low modes. This can be clearly seen from Figs. [13] 
and |14| This degeneracy indicates some higher symme- 
try that includes chiral SU(2)l x SU(2)ji as a subgroup. 
It is a highly exciting question what this higher symme- 
try is. It will be seen from the following subsection that 
baryons also point to some higher symmetry. 



[(0,i) + (i,0)] x [(0,i) + (i,0)] 
= (0,0) + (i,!) Q + (!,!)(, + (0,1) + (1,0) . (6) 
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FIG. 13. Summary plots: Baryon (l.h.s.) and meson (r.h.s.) masses as a function of the truncation level. 
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FIG. 14. Summary plots: Baryons (l.h.s.) and mesons (r.h.s.) in units of the p-mass at the corresponding truncation level. 



D. Baryon chiral multiplets 

If chiral symmetry is restored and baryons are still 
there they have to fall into (some of) the possible bary- 
onic parity-chiral multiplets. There are three different 
irreducible representations of SU{2)l x SU{2)r x d for 
baryons of any fixed spin: 

(i,0) + (0,i), (f,0) + (0,|), (|,1) + (1,|). (7) 

The first representation combines nucleons of positive 
and negative parity into a parity doublet. The second 
representation consists of both positive and negative par- 
ity A's of the same spin. Finally, the third representa- 
tion, that is a quartet, includes one nucleon and one Delta 
parity doublet with the same spin. 

Extraction of the chiral cigenmodes of the Dirac oper- 
ator leads to a systematic appearance of the parity dou- 
blets, as it is clearly seen from Figs. [13] and [Mj There 
are two degenerate nucleon parity doublets with the same 
mass. There are also two distinct A parity doublets, but 
with different mass. Since our interpolators have spin 
J = | for nucleons and J = f for Delta's, we cannot see 
possible quartets of the (a, 1) + (1, |) type. 

It is very interesting that the two nucleon parity dou- 
blets get degenerate, while the two Delta doublets are 
well split. The former hints at a higher symmetry for the 
J = I = | states, while this higher symmetry is absent 
for the J = I = | states. 



E. On the origin of the hyperfine splitting in QCD 

The A — N splitting is usually attributed to the hyper- 
fine spin-spin interaction between valence quarks. The 
realistic candidates for this interaction are the spin-spin 
color-magnetic interaction [SSI 121] and the flavor-spin 
interaction related to the spontaneous chiral symmetry 
breaking [55] . It is an old debated issue which one is 
really responsible for the hyperfine splittings in baryons. 
Our results suggest some answer to this question. Once 
chiral symmetry breaking is removed, which happens for 
the ground N and A states after extraction of the 50-60 
lowest eigenmodes, the A—N splitting is reduced roughly 
by the factor 2. With the restored chiral symmetry the 
effective flavor-spin quark-quark interaction is impossi- 
ble. The color-magnetic interaction is still there. This 
result suggests that in our real world the contribution of 
both these mechanisms to the A — N splitting is of equal 
importance. 



VI. CONCLUSIONS 

We have studied what happens with different mesons 
and baryons upon modifying the valence quark propa- 
gators by removing the lowest lying eigenmodes of the 
Dirac operator. These eigenmodes are directly related 
to the quark condensate of the vacuum via the Banks- 
Casher relation. Consequently, upon removal of the low- 
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est eigenmodes we artificially restore chiral symmetry, 
what we call "unbreaking" of the chiral symmetry. We 
study the evolution of the hadron masses with the num- 
ber of extracted lowest eigenmodes in dynamical lattice 
simulations. There are a few interesting observations. 

First, the quality of the signals from the hadrons af- 
ter removal of the chiral eigenmodes of the valence quark 
propagators become much better than with the untrun- 
cated propagators. Most probably this is related to the 
fact that we artificially remove the pion cloud of the 
hadrons. 

Second, from the spectral patterns both for the ground 
and excited mesons and baryons we conclude that con- 
finement is still there while the chiral symmetry is ar- 
tificially restored. Restoration of the chiral symmetry 
is evidenced by the fact that hadrons, both baryons and 
mesons, fall into different parity-chiral multiplets. At the 
same time there is a clear evidence that the broken U(\)a 
symmetry is not restored. 

Third, some distinct parity doublets get degenerate 
upon chiral symmetry restoration. This indicates that 
there is some higher symmetry in the chirally restored 



regime, that includes the chiral group as a subgroup. 

Finally, from the comparison of the hyperfine A — N 
splitting before and after unbreaking of the chiral sym- 
metry we conclude that in our real world both the color- 
magnetic and the flavor-spin interactions between valence 
quarks are of equal importance. 
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